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Abstract 

This work was mainly driven by the desire to explore, to what extent embedding some given geometry in a higher 
dimensional flat one is useful for understanding the causal structure of classical fields traveling in the former, in 
terms of that in the latter. We point out, in the 4-dimensional (4D) spatially flat Friedmann-Lemaitre-Robertson- 
Walker universe, that the causal structure of transverse-traceless (TT) gravitational waves can be elucidated by 
first reducing the problem to a 2D Minkowski wave equation with a time dependent potential, where the relevant 
Green’s function is pure tail ~ waves produced by a physical source propagate strictly within the null cone. By 
viewing this 2D world as embedded in a 4D one, the 2D Green’s function can also be seen to be sourced by a 
cylindrically symmetric scalar field in 3D. From both the 2D wave equation as well as the 3D scalar perspective, 
we recover the exact solution of the 4D graviton tail, for the case where the scale factor written in conformal 
time is a power law. There are no TT gravitational wave tails when the universe is radiation dominated because 
the background Ricci scalar is zero. In a matter dominated one, we estimate the amplitude of the tail to be 
suppressed relative to its null counterpart by both the ratio of the duration of the (isolated) source to the age of 
the universe rjoi and the ratio of the observer-source spatial distance (at the observer’s time) to the same rjo- In a 
universe driven primarily by a cosmological constant, the tail contribution to the background geometry a[?7]^77^^ 
after the source has ceased, is the conformal factor times a spacetime-constant symmetric matrix proportional 
to the spacetime volume integral of the TT part of the source’s stress-energy-momentum tensor. In other words, 
massless spin-2 gravitational waves exhibit a tail-induced memory effect in 4D de Sitter spacetime. 
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I. MOTIVATION AND INTRODUCTION 


The geometry of our universe appears to be well described by Einstein’s equations with a cosmological 
constant A = 


= SvrGNn. (1) 

At zeroth order and at very large scales, contains an isotropic and homogeneous background 

matter-energy distribution that drives the evolution of a 4-dimensional (4D) spatially flat Friedmann- 
Lemaitre-Robertson-Walker (FLRW) geometry parametrized by conformal time r] and 3 spatial coor¬ 
dinates X, i.e., 

diag[l,-l,-l,-l]. (2) 

At first order, + 5T^u also includes the inhomogeneous matter perturbations necessary 

to describe the finer structure present in the universe - clumping of dark matter, for instance - once 
cosmologists try to probe it at smaller scales and higher resolution. These perturbations will also 
produce inhomogeneities in the metric, so that now 

g^,u[v,x] = a[rf\^ -h h^t,['n,x]). (3) 

R is possible to perform a scalar-vector-tensor decomposition of both the matter and metric ha 0 
fluctuations, such that at linear order in these helds, eq. Q would yield separate partial differential 
equations (PDFs) for each of the perturbations transforming differently under the rotation group SO 3 , 
an isometry group of the background 4D FLRW geometry in eq. 

This paper is specifically about understanding the causal structure of the transverse-traceless (TT) 
metric perturbations 

h^ydx^dx'^ = Dijdx^dx^, 

with = 6 '‘^diDjk = 0. The Dij, describing gravitational radiation capable of carrying energy- 

momentum to infinity, is governed by 

UDij = levrGNn^J), (4) 

^ For a detailed and pedagogical treatment of cosmological perturbation theory, see for example, Chapter 5 of Weinberg’s 
cosmology textbook [T]. 
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with denoting the part of that is subject to the TT constraints, = 0.^ 


■jk 


Now, the □ in @ is the scalar 


one, i.e., 


UDij = 






i>Dij I , 


(5) 


where = a^, so that the contribution to the deviation from eq. Q due to the presence of II-J^ is 
given by 


Dij [r], x] = IOttCn 


dry' j d^xa[r]']'^G^[r],x-,r]',x']u‘^\r]',x‘ 


( 6 ) 


with the retarded G'^ obeying 


[??, X] r/', x'] = [r?, f; r/', f'] 


6 [r] — — x'] 

a[r]Y a[ri'Y 


(7) 


We see that the study of the causal structure of TT gravitational waves (GWs) propagating in our 
universe is the same as that of a minimally coupled massless scalar field. By causal structure, we are 
referring to the fact that, in a curved spacetime, particles that are otherwise massless in 4D Minkowski 
spacetime no longer travel strictly on the null cone - they travel both on and within it. (The part of 
the field traveling inside the light cone is usually called the “tail”.) While this is known in principle 
[2], a thorough understanding of the causal structure of waves propagating in physically important 
geometries is lacking, particularly in black hole spacetimes. The black hole case is not only a rich 
problem, because caustics can form from the focusing of null geodesics - indicating the causal structure 
of the Green’s functions in Kerr spacetime will likely be highly non-trivial - but understanding it has 
become fundamental to a successful prediction of GW signatures from Extreme-Mass-Ratio-Inspiral 
systems.^ 

In cosmology, photons and gravitons are assumed to follow null geodesics, and this is usually justified 
via JWKB arguments. If the universe is perfectly homogeneous and isotropic, this actually is an exact 
statement for electromagnetic fields because of the conformal symmetry enjoyed by the 4D Maxwell 
action. However, as we will see explicitly below, while light requires the inclusion of inhomogeneities in 
eq. Q to develop tails [8], GWs already do so in the 4D spatially flat FLRW background geometry of 
eq. (§ - except during radiation domination. Moreover, current and future generation GW detectors 
have the potential to listen in on sources at cosmological distances. One may wonder if the tail part 


These TT requirements provides 4 equations, reducing the 6 components of the symmetric Dij and of Ilh ’ down to 2 
independent ones. 

® See [3] through for a sample of recent work related to Green’s functions in Kerr/Schwarzschild spacetimes. 
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of such GW signals can thus acquire additional properties during their propagation, as a result of 
their interaction with the geometry of the expanding universe, such that they can in turn be used to 
tell us something about the history of cosmic evolution. The estimates made later on in the paper 
inform us that, unfortunately, this is unlikely, i.e., the dominant effect of cosmic evolution on Dij 
is that its waves propagate at unit speed and redshift as \/a[rj\. Nonetheless, we remark that the 
computation of GW signals from astrophysical sources, such as in-spiraling compact binary systems, 
are usually performed with asymptotically Minkowskian boundary conditions. A treatment of how 
GWs from astrophysical systems propagate over such cosmological distances would presumably have 
to involve drawing a correspondence between the stress-energy-momentum [r/, x] in eq. ([^ with its 
counterpart in the post-Minkowskian/Newtonian treatment, by demanding that the far field asymptotic 
predictions of the latter coincides with the near/intermediate zone ones of the former in equations eq. 

(§. 

One of the key technical goals of this paper is to continue the investigation initiated in [9] and [IHl , to 
see if embedding the relevant geometry in some appropriate higher dimensional Minkowski spacetime 
will aid in understanding the causal structure of waves in the former, since the Green’s functions 
encoding causal structure information is known explicitly in the latter in all dimensions greater or 
equal to 2. In the current cosmological context, we shall find it useful to exploit the spatial translation 
and rotation symmetries of the 4D FLRW geometry in eq. as well as a conformal re-scaling of the 
massless scalar/graviton Dij, to first reduce the problem to a 2D one. Following that, we will show how 
this allows us to cleanly separate the part of Dij that transmits signals on the light cone from its tail 
part. In particular, the light cone term of its Green’s function - the “direct” part, as it is usually known 
in the gravitational physics literature - is completely determined by the flat 4D Minkowski counterpart 
modulated by the conformal re-scaling 1 /[a\rf\a[r]']), where r] and r]' are the observation and emission 
time respectively. The tail part of Dij is what carries physical information about how the graviton 
interacts with the background spacetime. 

In this paper we will, for technical convenience, mostly deal with symmetric - retarded plus advanced 
- Green’s functions.^ Because the retarded/advanced Green’s functions are usually normalized so that 
they can be said to be sourced by some appropriately defined unit strength spacetime point source, 
the corresponding PDF for the symmetric Green’s function G[y, y'] is the field observed at y^ = ( 77 , x) 
sourced by a spacetime point source at y'^ = (t/', x') of amplitude 2. Specifically, the symmetric Green’s 

^ Retarded Green’s functions will be labeled with a ^ superscript, and is the symmetric counterpart multiplied by Q[r)—r)'\\ 
moreover, throughout this paper, 0[... ] is the Heaviside step function and (i[. •. ] is the Dirac delta-function. 
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function Gi\y,y'] equation for Dij obeys 


nyG^[y,y']=ny,G^[y,y'] = 2 


5[r] — — x'] 


a{r]Ya[ri'Y 


( 8 ) 


We will define mm 


, G4 

G4[x,x]^ , 

a[y]a[y'] 


(9) 


which will then lead us to the following 4D Minkowski wave equation with a time-dependent potential: 


with 


(52 - d^^d^dj + U[y]) G 4 = ( 52 , - + U[n']) G 4 = 25[ri - _ f'], 


rrr 1 _ ^[v] _ 

a[y]~ “ dr/2- 


( 10 ) 


( 11 ) 


For reference, the Ricci scalar TZ of the 4D spatially flat FLRW universe in eq. Q is proportional to 
this potential U: {a‘^/ 6 )TZ = —dja. 

In section (|TT|) , we will begin by explaining why any Green’s function equation of the form in eq. (10), 


with a space-independent but otherwise arbitrary potential - not necessarily arising from a cosmological 
context - can be reduced to a 2D one. We will then use this 2D —)■ 4D prescription to cleanly split the 
“direct” part of the cosmological Green’s function from its tail from the outset. In both 2D and 4D, 
the relevant wave equations and the light cone boundary conditions obeyed by the tail functions, will 
be discussed. We will apply what we have learned, to solve for the TT-GW Green’s function in a power 
law cosmology, and use it to estimate the size of the tail effect in a radiation, matter and cosmological 


constant dominated universe. In section (III), we offer a different perspective on this 2D wave equation 


with a potential, by showing how embedding this 2D world in a 4D flat one allows us to see that the 


2D Green’s function can be sourced by a 3D scalar field. We summarize our findings in section (IV) 
A significant portion of the analysis in sections 


and (III) is based on appendix There, we 

discuss the causal structure of the Green’s function associated with a flat spacetime 2D wave equation 
with an arbitrary potential. We will do so directly in 2D and also by embedding the 2D world in a 
flat 4D one. An infinite Born series solution is derived for the retarded 2D reduced Green’s function 
G 2 and for its tail function. The light cone boundary conditions obeyed by the first derivatives of this 
tail function are also worked out. In appendix (|^ we use position space methods to argue that the 
symmetric Green’s function of the 2D Minkowski wave operator d‘^ is given by 1/2 everywhere within 
the light cone of its source and zero outside of the null cone. 
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II. 4D CAUSAL STRUCTURE FROM 2D 


We will begin by explaining why any Green’s function G 4 in 4D that reflects the spatial translation 
and rotational symmetries of its background geometry, may be reduced to a 2D problem. Such a high 
degree of symmetry means both Green’s functions ought to depend on the spatial locations of the 
observers and sources through the Euclidean distances between them. Let us denote = a:, x'^ = x', 
x± = (0,x^,x^) and x'j_ = (0, x'^, x'^). Next, we proceed to consider 


G2[r],rj';E]= / d^x'^G^ 


r], 77 '; J(x - x')^ + (x_L - x'l)^ 


( 12 ) 


In what follows, it will be useful to work in terms of the symbol R. If it occurs within a 2D object, it is 
R = |x — x'l; if it occurs in a 4D object, then it is i? = |x — a?|. Similar remarks apply to Synge’s world 
function = (l/ 2 ')(y — y')^, which is half the square of the geodesic distance between the observer 
at y and the spacetime point source at y' in flat spacetime. By switching to cylindrical coordinates 
r_L = |al_L|, followed by p = eq. ( 12 ) then reads 


G2\n,r]']R]=2'K j dppGA[r},'n']p\. 

Jr 

Differentiating both sides with respect to R then tells us G 4 does indeed follow from G 2 : 

1 dG2[rj,v'-,R] 


(13) 


G4[v,v-,R] = - 


2ttR 


dR 


(14) 


^The 2D Green’s function, in turn, obeys the 2D analog of eq. (10), 


-dl + U[p\) G 2 = - dl + U[ini']) G 2 = 25[r] - ini']5[x - x']. 


(15) 


To see this, we apply W ’2 = + U[r]] to both sides of eq. ( 12 ). Interchanging the differentiation 

and integration on the right hand side, and then adding and subtracting the Euclidean Laplacian 


= 02+^3 with respect to x±, and finally invoking eq. (10) obeyed by G 4 , we have 


[d'^ - dl + U[r]]) G 2 [i],r]'; R] = 2S[r] - p']d[x - x'] + d^x_LV^G4 p, p'; J {x - x'Y + 

/ Tn>2 * 


• (16) 


We recognize the remaining integral on the right hand side of eq. (16) to be a surface term at spatial 


infinity. As long as the G 4 respects causality, then this surface term must be zero because the source 


® This result when applied to Minkowski Green’s functions, is why there exists a recursion relation allowing one to 
construct all even dimensional ones from the 2D case and the odd ones from the 3D result; see m for a pedagogical 
discussion. 
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and observer, for fixed times and x, x', must lie outside each other’s light cone as |x_l| —)• oo. (Similar 
manipulations hold for the equation with respect to (ry',x').) 

Because of the primary role played by the 2D wave equation (15), we will in appendix study its 
properties for a general potential = (??, x)] - the reader is invited to take a detour there before 

returning to the main body of this paper. Here, we will summarize the main results, but specialize to 


the case where U depends only on the conformal time rj. The solution to eq. (15) takes the general 
pure tail form 


G 2 [r?,x;? 7 ',x'] = ^ 0 [cr] J[? 7 , r/', i?]. 


(17) 


where J itself obeys the homogeneous 2D wave equation 

(52 _q 2^ j ^ ^ j ^ 0^ (Ig) 

and the light cone boundary condition J[a = 0] = 1. These facts, together with —{l/R)dR = d/da and 


eq. (14) immediately imply the 4D reduced Green’s function itself takes the form 

04 = G4M + , 

where G 4 [fj] is the 4D minimally coupled massless scalar symmetric Green’s function, 

5[d] 


G4[d] = 


dvr 


(19) 


( 20 ) 


Without solving any PDEs, we have managed to isolate the light cone part of G 4 from its tail. Note 
that, because 4D FLRW cosmology is conformally flat, its light cone is the Minkowskian one 


a = 0. Witness too that the 4D tail term dJ/da in eq. (19) is entirely determined by the homogeneous 
solution to the 2D wave equation. 

Because it may be useful for 4D cosmology, we record here that the analogous homogeneous equations 
for the 4D graviton tail in terms of rj, r/' and a treated as independent variables. 


dr^ {2{r} - ■q')daJ') + d^ {2adaJ') + (5^ + U[rj\)J' = 0, 
dn' ( 2 ( 77 ' - r])daJ') + da (2adaJ') + (5^/ + U[iq'])J' = 0. 


( 21 ) 


In eq. (21) and the following paragraphs, J = and J = daJ- As one may expect, evaluating 
a = (l/ 2 )((r 7 — 77 )'^ — (x — x')^ — (x_l — x(j_)^), we find that eq. ([^ is simply the homogeneous version 


( 22 ) 


of the 4D wave equation for the reduced Green’s function (eq. ®), 

(92 - S^^didj + t/[r 7 ]) J' = (4 - d^^di^djf + U[ri']) J' = 0. 
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The light cone boundary condition for J' with an arbitrary but space-independent potential U[r]] is 

/ dJ[r],ri';a = 0]\ 1 




da 




2 {r] - i) 


dfU[f]. 


(23) 


In a 4D spatially flat cosmology, U = —d/a and hence eq. (23) yields 

'dJ[r],r]';d = 0] 


da 


drj' 




r^,n' 2 (r? - r]') 

In terms of W 2 J = {d/j — dj^ + U) J, eq. (21) follows from eq. (A3) and 


a[r ]"]' 


(24) 


0 = d^W 2 J = (W 2 + 2d^)J'. 


(25) 


As for the boundary conditions in eq. (23), it is convenient to use the light cone coordinates = rjdzR 


and = r]', so that rj = (^+ + ^ )/2 and a = (1/2)(^+ — r/')(C “ V)- This in turns implies 

d±J=lj+l{^^-vV'. 


(26) 


We now consider the light cone limit u = 0 by setting —)• 77 '. Because J[r], r]';d = 0] = 1 for all times, 
we must have J[d' = 0] = 0. Application of equations (A26) and ( A27| ) to eq. (26), and recognizing 
v" = (?^^^ + returns 

1 r(?^+»7')/2 

/ T}' 


^- - - r}'jj'[q,r}'-,a = d] = -- J ^ dri”U[ri”]. 


(27) 


The are arbitrary at this point, but we would like results expressed solely in terms of r] and 77 ', so 
we may put + 77 ' —>■ 2r] and thereby verify eq. (23). 


What eq. (19) and eq. ([^ inserted to eq. (|^ teaches us, is that a physical source in any spatially 
flat FLRW universe produces a null GW front that takes a universal form: 

,nSP [r]r,x'] 


T>-7^[77,f] = 4 Gn / d^f'a[r7,.] 


i3 _ 

a[ 77 l|x — x' 


(28) 


where rjr ^ rj — \x — x'\ is the retarded time. If the source is isolated, we may define if = (l to be in its 
interior. Then in the far field limit - specifically, whenever |i:| is much larger than the spatial extent 
of the source - we may deduce this null portion of the GW detected at ( 77 , x) is, roughly speaking, the 
spatial-total of the stress-momentum of the source at retarded time, modulated by the inverse physical 
distance between observer and source: 




a 77 X 


d^x'a [77 — |i;|]^ II^^J^ [77 — |x|, if] . 


(29) 






















FIG. 1: Causal structure of TT GWs. This spacetime diagram depicts a hypothetical astrophysical process, 
where a massive star undergoes core collapse and goes supernova (right world line). The dashed-dotted segment 
of the right world line represents the full duration during which GWs are produced: before that, the collapse has 
not started; after that, the system has settled down completely. These GWs are detected by a distant detector 
(left world line). We assume that the background geometry is that of a 4D spatially flat FLRW universe. The 
black dashed lines emanating from the detector’s world line are the past light cones of events B and C. 
The bottom pair of light gray dashed lines emanating from the right world line is the forward light cone of 
the beginning of the collapse; the top pair is that of the end of the process. The light gray shaded region of 
spacetime is filled with GWs propagating both on and inside the null cone. The darker-gray region of spacetime 
is filled with GW tails only. The detector at A sees no signal. The detector at B sees a “direct” signal from B' 
but also a tail from the dashed-dotted segment of the massive star before B'. At G the signal received is the 
accumulation of the GW tails from the entire dashed-dotted segment. As we show in the main text, in a radiation 
dominated universe, there are in fact no tails, so the detector sees nothing at C. In a matter dominated universe, 
the signal at C is independent of the spatial location but decays with time in an expanding universe. In a de 
Sitter background, the Green’s function tail is a constant; the detector that was operational from A through C 
would sense a permanent change in Dij that is proportional to (H/Mpi)'^ and to the total contained in the 
dashed-dotted segment, i.e., (H/Mpi)"^ J (f^x'^y\g[r]']\ 11 [ 77 ' ,x']. 

Unlike its direct part, the knowledge of tail dJjda requires solving its full wave equation. Hence, it is 
the part of the GW moving slower than unit speed that actually encodes information regarding how 
these massless spin-2 particles interact with the background geometry of eq. Q. 
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We now turn to solving the GW tail in a power law cosmology. 

Power Law Cosmology A scale factor that is a power law in conformal time, 

V 


a[v] = ^ 

\Vo 


( 30 ) 


can be used to approximate the major epochs in cosmology - inflation/dark energy/de Sitter (p = — 1), 
radiation (p = 1) and matter domination (p = 2). This means U = —aja becomes 

p{p - 1) 


U[v] = -- 




(31) 


and the corresponding 2D wave equation for J in eq. (18) can be reduced to an ordinary differential 
equation. (We would like to attribute this realization to Nariai [13 •^) Let J depend on spacetime 


solely through the combination s = aj{rir]'). Then eq. (18) becomes 


^ {s(s + 2)J"[s] + 2(s + l)J'[s] - {p- l)pJ[s]} = 0. 

The regular solution obeying the light cone boundary condition J[a = s = 0] = 1 reads 


J[s] — P-p[s + 1] — Pp_i[s + 1] — 2 F 1 


1 -p,p;l;--j 


a 


s = 


rjT]' 


(32) 


(33) 


We have provided, for the reader’s reference, 2 equivalent ways of expressing the Legendre function P^, 
as well as its hypergeometric function 2 P 1 representation; see the p = 0 limit of equations 8.702 and 
8.704 of |15j . To sum, the symmetric Green’s function G 2 obeying 


.2 pip - 1 ) 


q2 _ q2 _ ^) \ P< ^ __ „/ _ / o2 a2 pip 1) \ Pi r__ 


7]^ 


G 2 [ri,x-,r] ,x] = [d , - d^, - 


rj' 


n 


G2[ri,x-,r],x 


= 25[r] — ri']6[x — x'], 


(34) 


IS 


G2[?7 ,x;??,x'] = -Q[a\P-p 


a 


a 

^ -7 

rjTjf 

= 



(35) 


with = {(p — rj'Y — ix — x')^| /2. At this point, putting back the conformal re-scaling (eq. ([^), we 
have (re-)derived the 4D symmetric Green’s function of the minimally coupled massless scalar - and 
therefore, also that of the transverse-traceless graviton (obeying eq. (j^) - in a spatially flat FLRW 


background power law universe (eq. (30)): 

1 


Gi[n,v',R] = 


ATra[r]]a[r]'] 

1 

47ra[7?]a[r/'] 


(5[cr] -I- 0[cr] 


P'-p [l + (T/{r]ri')]\ 


rjT]' 


( 36 ) 




2 -p,p+l;2;- 


a 


2r/r/' 


I also wish to thank Dai De-Chang for bringing m to my attention. 
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where the prime on the Legendre function denotes derivative with respect to its argument, and here 
a = (l/2)((r/ — rj'Y — {x — For ease of comparison to Nariai’s results in eq. (3.5) of [13], we also 

display the solution in terms 2 .^ 1 .^ One can check that, since 2 F"i[a, /3; 7 ; 0] = 1, the tail portion of the 


solution in eq. (36) satisfies the light cone boundary condition in eq. 

5J[7, 7 '; fj = 0 ] p{p — 1) 


da 


2rjr]' 


(37) 


Radiation domination The second line of eq. (36) makes manifest the fact that there are no GW 
tails in a radiation dominated (p = 1) universe. Mathematically this can be understood from U in eq. 


(31) being zero, and therefore G 2 = 0[ct]/2. One may also recall that, if the radiation were primarily 
photons, because the Maxwell action is conformally invariant in 4D, the trace of its stress-energy- 
momentum tensor is zero - and so is the background Ricci scalar TZ, by Einstein’s equations. That U 
is zero then follows, as we have previously observed that it is proportional to the Ricci scalar. 

6[a] 


G 


(Radiation) 


[ri,V-,R] = 


'n 

. r 1 r /i’ ^ (0,oo), a[rj] = —. 

47 ra[r]\a[r]'\ 770 


(38) 


(T) 

Therefore the contribution to the geometry Dij given the TT part of some source IT) •' is completely 

J tj 


determined by the universal form in eq. (28). 


^ We mention here that Caldwell in m had claimed to have solved the Green’s function of the TT graviton in a FLRW 
universe with arbitrary spatial curvature. We do not believe his results are correct. When K = 0 (spatially flat FLRW), 
the V/{ri — rj') in Table I of [16] corresponds to our tail dJjda. However, on the very first line of Table I, the power 
law cosmology TT-GW tail obtained there was dJ/da = p{p — l)/{2pr]') - Caldwell’s is our (p,rj') - and is 

therefore inconsistent with the solution obtained by m, [m and in eq. (|36||. Indeed, his spatially flat tail result is 


really the light cone boundary condition in eq. (37l; compare his (III.6) or (III.7) with eq. (24l. One potential source 
of error is that Caldwell appears to have incorrectly asserted that Synge’s world function in FLRW spacetime is given 
by O' = (1/2)(R^ — iv ~ v')^) ~ S66 the statement right before his eq. (III.6) - which he then used in his eq. (11.7) to 
calculate the TT GW Green’s function. 

For the reader’s reference, we compare the methods used here, m, and m, to obtain the minimally coupled massless 
scalar Green’s function in a 4D spatially flat power law cosmology. In eq. (2.13) of [TH] and eq. (3.8) of [TT], a 
Fourier transform/mode-integral representation involving Bessel functions was found for the Green’s function. Nariai 
m replaced the Bessel functions with their asymptotic expansions and proceeded to work out a power series expansion 
for the Green’s function tail in eq. (2.17). Following that, an ansatz for the exact solution was proposed in eq. (3.1), 
which lead to a hypergeometric (ordinary) differential equation in terms of our al{2r\r{'). The appropriate solution of the 
two linearly independent ones was then selected by matching it onto the power series in eq. (2.17). In Haas and Poisson 
m [S] was consulted to directly evaluate the mode integrals in terms of the Appell function 77; see their Green’s 
function tail in eq. (3.10). Hence, their solution for a generic power law cosmology is more complicated that Nariai’s 
m and ours (eq. ( |36| )) in that - by referring to [15], for instance - F 4 appears to have fewer tabulated properties than 
2 F 1 . On the other hand, in this paper, we provided a largely self-contained derivation without resorting to the Fourier 
transform. We showed that the 4D symmetric Green’s function is simply related to the first derivative (with respect 
to a) of its 2D cousin. The 2D Green’s function, in turn, not only takes a pure tail form, the tail J is fixed uniquely 
because it satisfies the homogeneous 2D wave equation and is unity on the light cone. To be sure, however, note that 
our 2D power law solution was gotten by postulating that the tail J was a function solely of a/{rip'), and this was 
certainly predicated on Nariai’s results. 

Note that Nariai uses observer time t, where a\t\ oc F and his background FLRW geometry is AR — a[t]'^dx‘^', our p is 
related to q there via p — q/{l — q); Nariai’s v = (3g — 1)/(2(1 — q)), t = p — p', ^ = 2a and « = al{2pp'). 
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By setting p = 2 and p = — 1 in eq. (36), one can verify that the reduced 4D Green’s function 


G 4 of the TT graviton is in fact the same object in both matter and cosmological constant dominated 
universe. Their difference in G 4 is entirely due to the overall conformal re-scaling 1 /{a[r]]a[ri']). 

Matter domination The matter dominated TT graviton Green’s function tail is constant in 

3 3 

space but goes as l/(a[ 7 ?] 2 a[r/'] 2 ), 

2 




1 




4 

/I I r 1 r n T © [d'] 

47 r \a[ri\a[r] \ 

and the deviation from eq. © due to ,x'] is given by 


p,p' G (0,oo), a[p] = ( ^ ) , (39) 


Di^ = 

U U ’ 


(40) 


where the wavefront that travels at unit speed given by eq. (28) and the tail part of the GW is 

fr;—|x—a;'|—0+ 


D 


(tail) 


[V, x] = 


4Gn 

2 

7?2a[p]2 


f 


d77'a[r7']2 / [^/^^] • 


(41) 


Suppose the strength of the source peaks at time 77 *, and suppose At r\j /peak width dVa[r?'] is the physical 
duration of the source in the cosmic rest frame. The direct part of the signal in eq. ( |28[ ), in the far 
field limit, is then roughly bounded by 


dI]\p,x] 


^ 4Gn 

a[77]|x| 


whereas the tail is roughly bounded by 




(tail) 


[q,x\ 


4Gn At 


pla[pya[p^ 


d^f'a[77*]^n^P [p^,x' 




(42) 


(43) 


Suppose we normalize the scale factor such that r/ = po is the observer’s time; up to factors of unity, and 
assuming matter domination throughout cosmic history, r/o is also the age of the universe. According 
to the observer, the ratio of the peak tail amplitude to that of the null cone signal can be estimated as 


MtaiX) f 


D^^\pq,x\ 




\m ) \m 


(44) 


The tail effect, in a matter dominated spatially flat FLRW universe, is suppressed relative to its null 
cone counterpart, by the ratio of the duration of the source to the age of the universe times the ratio 
of the observer-source distance (at the observer’s time) to the age of the universe. The way to beat 
this suppression is to have the source of GW reside at a very early epoch of this universe, so that 
a[v*] ^ i{^t/po){\x\/po)f/^. 
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de Sitter 


The de Sitter TT graviton Green’s function tail is a constant: 




(5 [a] 


+ 


0[d] 


G (-00,0), a[rj] = ^. 


(45) 


47r \^a[7y]a[r7'] ' rj^ 

(This is consistent with the de Sitter minimally coupled massless scalar result in [9].|in].) Therefore 
the metric fluctuation Dij = engendered by the TT part of the source has a 


null cone portion given by eq. (28) and a tail part that reads 


T), 


(tail) 


[V, x] = 


4Gn 


frj—|x—x'|— 0 + 


dr]' f d^x [n',x' 


(46) 


Vo J—ao 

/rpN 

Assuming lib [r]',x'] describes a source that radiates GWs over a finite duration, any detector that 

was present before, during, and after the GW train has passed the observer’s location will find that 

Dij does not decay back to zero but instead suffers a permanent “DG” shift AiDjj proportional to the 

(T) 

spacetime volume integral of lib >; 


AD - ^ 

^■^ij — 2 

Vo 


r0~ 




dr]' f d^fVl5[f/']|n-p [v',x' 


(47) 

(48) 


®For de Sitter spacetime, we may identify r]o = —1/H to be the negative reciprocal of the Hubble 
expansion parameter. This means the size of this memory effect is governed by the square of the 
ratio of Hubble to that of the Planck mass (77/Mpi)^, with Gn ~ l/Mpp while the magnitude of 
the components of the constant symmetric matrix Hij is, heuristically speaking, the TT-part of the 
spacetime-total of the stress-momentum of the source. 

Just as for the matter dominated case, let us compare the peak amplitudes of the direct signal to 
the tail. Again, we assume H^J^ itself peaks at ry*, and the source duration is At rsj /peak width [?/']■ 
Then - recalling “ we may estimate that in the far zone 


D 


dI]\v,x] 


^ 4Gn 

^ a[v]\x\ 


4Gn At 
Vo 


f d^x'a[r]^fll^^^ [v*:^] 

Jr3 

[ d^f'a[ry*]^n^^p [7y*,f'] 

JR3 


(49) 

(50) 


® For foundational work on the gravitational memory effect in asymptotically flat spacetimes, see m and [H]; for more 
recent theoretical investigations, [HHS]; for observational searches, see for instance [^. In addition, recent work in 
|23| and have connected this asymptotically Minkowskian memory effect to the low frequency limit of the Ward- 
Takahashi identities obeyed by graviton scattering amplitudes due to the Bondi-van der Burg-Metzner-Sachs (BMS) 
symmetry at null infinity, commonly known as Weinberg’s soft graviton theorem. We pose the analogous question here: 
is the memory effect in eq. (471 related to some symmetries associated with asymptotically de Sitter spacetimes? (Note 
that de Sitter asymptotics are not as well studied as flat ones - see [25].) 
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The ratio of the tail amplitude to that of the null cone signal is roughly, [H ■ At){H ■ a[ 7 /]|x|), where 
a\x\ is the co-moving source-observer distance evaluated at the observer’s time. 

We illustrate what we have found regarding the causal structure of FLRW TT gravitons in Fig. Q. 
We reiterate that the estimates performed here provide strong evidence that, for the most part, TT 
GWs in our universe propagate on the null cone. The predominant effect of cosmic evolution on them 


is the redshift l/a[ry] in eq. (28) 


III. 3D SCALAR PERSPECTIVE: 2D EMBEDDED IN 4D 


We will now take a different perspective on the 2D wave equation of eq. (15), by demonstrating 


that it is possible to view the 2D world as embedded in some 4D Minkowski spacetime, such that the 
G 2 is sourced by an appropriate 2D time dependent plane source J[ri,ri']r\ = x^]: 

1 


G 2 [r],r]'-R]= j d?xxGi[a]J[r],'q'■x\], <y = 7 . {{'H - - x]) 


(51) 


where we have let the spatial dimension of the 2D world pierce the 2D x_L-plane at its origin x_l = 0_l, 


and (74 (from eq. ( 20 )) itself obeys 


dlGi[a] = 5y,G4[cr] = - y']. 

As it turns out, J can be viewed as a cylindrically symmetric scalar field living in 3D, obeying 

(52 + Vi + U[y]) J + (r? - r/')Vx ln[ri] • Vxd^J = 0, 

(d^, + Vi + U[y']) J + {y' - y)Vx ln[ri] • Vxdr,'J = 0. 

and the boundary condition 


(52) 


(53) 


= ()±] = 1 , 


V??, rj'. 


(54) 


We discuss in appendix the more general case where the U depends on both rj and the first spatial 


coordinate x; we refer the reader to it for the derivation of equations (53) and (54). Notice from the 


prescription in eq. (51), that since G 4 itself is Poincare invariant, we have attributed to the source J 


all the time-translation symmetry breaking effects encoded in G 2 , due to the presence of the potential 


U[r]]. The representation in eq. (51) also illuminates the causal structure of a signal due to a spacetime 
point source at ( 77 ', x') in the 2D world, in terms of that of the collective 4D signal due to the 2D plane 
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FIG. 2: ID space embedded in 3D space. 


This figure illustrates the causal structure encapsulated in the 


integral representation of G 2 in eq. (51). The large shaded oval is to be viewed as a portion of the infinite 2D 
plane source J[r],ri 'which comes into existence for only an instant at 77 '; whereas the observer’s time is rj, 
i.e., the elapsed time between observation and emission is rj — rj'. The ID space of the 2D world whose scalar 
waves are described by G 2 , pierces the 2D plane source J orthogonally at its origin (denoted by O). The 2D 
observer is located R away from O along the ID line. (In 2D R = \x — x'\, where x and x' are, respectively, 
the spatial locations of the observer and the source of the Green’s function. In 4D R= \x — x'\, with similar 
meanings for x and x' .) Because the ambient 4D spacetime is Minkowskian (the G 4 in eq. ([^), for a fixed R 
and elapsed time \r] — rj'l, no signal from any part of the 2D plane source J can reach the observer whenever 
\r] — r]'\ < R. Suppose instead {r] — r]'\ > R, then because massless scalars in 4D Minkowski travel strictly on the 
light cone, the observed signal - from this 4D perspective - receives contributions only from the (dotted) circle 
on J of infinitesimal thickness drj_ and radius defined by \ri — = R^ + r\. But from the 2D point of view 


\rj — ri'\ > i? is simply the statement that the observed signal is the tail of G 2 ; from eq. (55) we also see the 2D 
signal has no (5-function impulse at [77 — 77 '! = R. Finally, since it is the integrated signal that is being observed, 
there is no need to allow for the source J to have an azimuthal dependence. 


source J[r 7 , rj'; - we explain it in some detail through Fig.(|^.^^ Furthermore, by using eq. (20) and 


It it worth highlighting that the embedding perspective here should not be taken too literally because the source 
J[ 77 , 77 '; rj.], although turned on only for an instant, somehow needs to take into account the time of observation rj to 
properly fool the 2D observer she is detecting the field of a spacetime point source at x' . This “acausality” of the 
source never arises when solving Green’s functions within a given geometry; for example, in \I1G[X,X'] = &^'^^\X — 
X'\/ ^| 5 [X]g'[X']| the source location X' is completely independent of the observer’s X. 
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switching to cylindrical coordinates in eq. (51), we may demonstrate explicitly that G 2 [ri, r(\ i?] is pure 


tail (see eq. (A15)): 


G 2 [r], f/'; R] = r]'; rl = 2a]. 


(55) 


Utilizing —{l/R)dR = da and equations (14) and (20) on eq. (55), we again recover the null-cone versus 


tail split in eq. (@. Moreover, when the relationship r_i_ -H- 2a is kept in mind, a direct calculation 


will show that equations (53) and (54) is in fact just a re-phrasing of the 2D homogeneous equation 


obeyed by J[r/, a], namely 


[r,, ri-R\ - ^ + [/[??] J[7y, r?'; R\ 

= _ j(o,o, 2 ) ^ U[r^']J% r?'; i?] = 0 


(56) 


and the light cone boundary condition J[ri, rj'-, = 2cj = 0] = 1 - except here R = \J{rj ~ v'Y ~ ^1- 


Is there some general expansion scheme that would allow us to solve eq. (53) using, for e.g., a 


separation-of-variables technique, while incorporating the boundary condition eq. (54)? At the moment 


we merely pose the question, but will show how the power law universe results of the past section follows 


readily from equations (53) and (54). 


Power law cosmology From the 3D scalar point of view, when the scale factor of the universe 
is a power law in conformal time a[r]] = {tj/tjoY, we use the ansatz J[s = r^/(2ryr/')], and will find J 


obeys, via eq. (53), 


/ 2„2 




(s(s + 2)j"[s] + 2{s -I- 1) J'[s] - {p - l)pJ[s]) = 0. 


(57) 


This is, of course, equivalent to eq. (32). Requiring J[r_L = 0] = 1, eq. (54), leads to the solution 


J[r],'n',r±] = P-p[s -h 1] = Pp-i[s -1- 1] = 2 F 1 


1 -p,p;l;- 


2J 


s = 


2 r]r]' 


(58) 


We recover the results of the previous section once we invoke eq. (55). 


IV. SUMMARY, DISCUSSION, AND A 5D QUESTION 

In this paper we have shown that, after the (known) conformal re-scaling of the TT graviton Green’s 
function in a 4D spatially flat FLRW universe, an understanding of its causal structure can be achieved 
by first reducing the problem to a 2D one. It turns out that the reduced 4D Green’s function’s “direct” 
part is equal to its Minkowski cousin, for any background cosmic history a[r]]. It is really the tail part 
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of the TT graviton that has physical information about how these massless spin 2 particles interact 
with the background FLRW geometry, since it is here that the wave equation needs to be solved in 
detail before the tail can be known explicitly. What we have uncovered is the 4D graviton tail function 
can be thought of as the derivative of the tail of the associated 2D Green’s function with respect to 
a, Synge’s world function in Minkowski spacetime. It can also be viewed as the radial derivative of 
a cylindrically symmetric scalar field residing in 3D. This 3D perspective comes from embedding the 
2D world in a 4D Minkowski spacetime. The 2D Green’s function is sourced by some instantaneous 
(at rj') 2D plane source J[rj,r]'] |x_l|], which intersects the ID space of the 2D world orthogonally at 
x'. This J is a 3D scalar because it depends on both time and the 2D space parametrized by x_l. 
Moreover, while in the 2D context the tail function obeys boundary condition on the light cone - it is 
unity there - the analogous boundary condition for the 3D scalar J translates to a spatial one, where 
|fj_| = 0] = 1. 

We have built on earlier work by Nariai m, and used our 2D and 3D pictures to recover the exact 
TT graviton Green’s function in a power law cosmology. We found there are no TT GW tails in a 
radiation dominated universe. In a matter dominated cosmology, the TT GW tail is constant in space 
but decays with increasing time (in conformal coordinates). We estimated, in this case, that the ratio 
of the peak amplitude of the tail part of the signal generated by an isolated source, to that of the 
direct part, to be proportional to the ratio of the duration of the source to the age of the universe 
and also to the observer-source distance (at the observer’s time) to the same. In de Sitter spacetime, 
the tail of the TT graviton Green’s function is a constant, and this leads to an inside-the-null-cone 
memory effect proportional to {H/Mpi)'^ ~ and to the TT-portion of the spacetime-total stress- 

momentum of the source. On the other hand, in 4D Minkowski spacetime, the massless scalar and the 
spin—2 graviton travel strictly on the light cone, i.e., they have no tails. As far as causal structure is 
concerned, therefore, the de Sitter memory effect uncovered here has no analog to its counterpart(s) in 
4D asymptotically flat spacetimes. 

As part of our explorations, we have studied in some detail - see appendix © - the properties 
of the 2D Minkowski wave equation with an arbitrary potential. This includes not only the pure tail 
nature of the corresponding 2D Green’s function; the light cone boundary conditions obeyed by the 
tail function J and its first derivatives; an infinite Born series solution for both the retarded version 
of the 2D Green’s function and for J; as well as how the 2D tail J can also be viewed as a 3D scalar 
held by embedding the 2D world in 4D Minkowski. As a by product, we have found that there are 
at least 3 other ways of obtaining the 2D Minkowski massive scalar Green’s function, apart from its 
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Fourier integral representation: by reducing the associated homogeneous wave equation into an ordinary 
differential equation; by switching to the 3D scalar picture, where the tail obeys a Laplace-Helmholtz 
equation; and by evaluating the tail’s infinite Born series expansion term-by-term. 

Dimension reduction, together with embedding, has allowed us to write down a PDF solely for the 
massless scalar/TT-graviton tail, before it is solved for explicitly. This is to be contrasted with what 
is usually done in cosmology, where a Fourier expansion of Einstein’s equations linearized about the 
background geometry of eq. ([^ is performed from the outset. Fourier space is crucial, of course, for 
applications such as predicting the power spectra of metric fluctuations that could be traced back to 
the inflationary epoch of the early universe. But to obtain information regarding the causal structure 
of TT-GWs, one would have to inverse Fourier transform the mode expansion of the Green’s function, 
and except for special cases where the mode functions are known explicitly - and simple enough for the 
integrals to be performed momentum space actually appears to hide both the general features of the 
light cone versus tail separation and the detailed solution of the tail itself. In fact our recognition here, 
that the Green’s function of the massless scalar field in 4D spatially flat FLRW can be cleanly split into 
its direct and tail pieces, can also be found in past literature - see for example, eq. (3.3) of [H] (which 
was a follow up of [I2j)- However, in where the case of a power law cosmology a[r]] = (r//?7o)“ 
was solved, the tail part of the Green’s function was still gotten via the Fourier transform of the whole 
Green’s function (equations (3.4), (3.6) and (3.8)). We see that working in position spacetime as far as 
possible, and utilizing the dimension reduction/embedding perspective, have let us arrive more readily 
at a simpler form of the solution consistent with that found by Nariai [13] (See footnote ([^ for the 
more subtle details.) 

We actually started out this 4D spatially flat FLRW investigation motivated by the fact that our 
4D universe can be embedded in 5D Minkowski. Is there a “line mass/charge” in 5D piercing the 4D 
FLRW world at x', such that the 4D observer is fooled into thinking that she is detecting the field 
generated by a spacetime point source at x', i.e., that of a Green’s function Gi[x,x'\I This is still 
unanswered. We hope, however, that the insights we have gained in this paper would help us further 
this goal. 

It could be that some integral transform or expansion scheme is useful for the solution of the tail part of the cosmological 
scalar Green’s function, but this we advocate, should be done directly for the PDE obeyed by the 2D tail itself, not the 
entire 4D Green’s function. Also, since the tail portion of the cosmological Green’s function, when reduced to the 2D 
context, is unity on the light cone, we suggest here that any useful method should be able to incorporate this boundary 
condition in a transparent manner. 


18 



V. ACKNOWLEDGMENTS 


I had exchanges with numerous people on the GW tail effect in cosmology. A partial list includes: 
Thorsten Battefeld, Lasha Berezhiani, Chiara Caprini, Dai De-Chang, Eanna Flanagan, Abraham 
Harte, Mudit Jain, Justin Khoury, Eric Poisson, and Vitaly Vanchurin. Much of the analytic work was 
carried out with Mathematica |2^. This work began while I was a postdoctoral fellow at the University 
of Pennsylvania, supported by NSF PHY-1145525 and funds from the University of Pennsylvania. Much 
of the analytic work performed in this paper was carried out using Mathematica |26| . 


Appendix A: 2D Minkowski Wave Equation With An Arbitrary Potential 

In this section, we wish to examine the dynamics of a scalar field obeying a 2D Minkowski wave 
equation with a potential U which is both time- and space-dependent. The dehning equation for the 
symmetric (retarded plus advanced) Green’s function (G 2 [C,C'] = G 2 [^^^]) reads 


(9| + U[e]) G2[?,?'] = (4 + [/[?']) G2[C,e'] = 


(Al) 


^^2 = {r],x) and = {rj',x') denote coordinates in (1 -|- 1)D; the subscript 2 in ^2 indicate the 
spacetime indices run from 0 to 1, so 9^ = — ^ri ~ 1 = d^, — d'i,, for instance. The 


key observation is that the solution to eq. (Al) takes the generic pure tail form 


where J obeys the homogeneous wave equation 




(A2) 


(4 + c/[e]) = (4 + me]) mee] = o, 


(A3) 


and the boundary condition that it is unity on the light cone 


J[cr = 0] = 1. 


(A4) 


Here, a is half the square of the geodesic distance between ^ and e 2D Minkowski. Because we are 


dealing with the symmetric Green’s function, The step function in eq. (A2) tells 


us scalar waves in 2D obeying eq. (Al) travel strictly inside the cone of its physical sources. This is 


to be contrasted with the 4D case, where there is an additional term proportional to J[d'], telling us 
scalar waves produced by a spacetime point source also receive contributions from a impulsive shock 
wavefront traveling exactly with unit speed. 
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We will now proceed to derive equations (A2) through (A4) by a direct calculation. Insert eq. (A2) 


into eq. (Al). On the left-hand-side, finds 


{d^ + U)G2 = ^ (d^O ■ J + 2d[a]d^ad'^J + O ■ (d^ + U)J) . 


We will argue in appendix (|B|) that, in Cartesian coordinates. 


(A5) 


a20[a] =40[u] 


(A 6 ) 


Thus, the first term from the left of eq. (A5) returns the desired d-functions in the right hand side of 


eq. (Al) if J = 1 when ^ Moreover, for eq. to hold, the coefficients of the <5- and 0 -function 


in eq. (A5) needs to vanish. Noting d^a = (^ — the d-function coefficient informs us J must be 
constant on the light cone: 


(e - i'rd^J[a = 0] = 0. 


(A7) 


But since J is unity at the apex of the light cone ^ that means it must in fact be unity everywhere 


(T = 0, i.e., we have eq. (A4). For the 0 function term to vanish everywhere inside the light cone, we 


see J must obey the homogeneous equation in eq. (A3). 

Massive Scalar in 2D Minkowski As an application, let us derive the massive scalar Green’s 
function in 2D Minkowski, 


(d^ + u) G 2 - e'] = 2<5(2) [e - , U = m^. (A8) 

We will assume, because of the highly symmetric nature of the problem, that J = J[my/^] depends 


on ^ and solely through a alone. The homogeneous equation (A3) then reads 


J"[x] + ^ + J[x]) = 0, 


y = my la 


(A9) 


The solution that obeys the light cone boundary condition J[x, u —>• 0] = 1 is the Bessel function 
jQ[m \/— ^0^]- We have, therefore, the solution to the massive scalar symmetric Green’s function: 


G2[^-e'] = ^e[d]Jo 


my la 


(AlO) 


This result can be cross-checked by performing the Fourier integral 

e[ri-v]G2[^-^']= ,2 . 2 ^ 2 ^ 

Jret 

where the retarded contour needs to be chosen for the /co-integral. 


(All) 
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2D Embedded in 4D Let us imagine that our 2D world is embedded in 4D Minkowski, and 
ask if there is some source in 4D that could in fact yield G 2 ■ Employing the minimally coupled massless 


scalar 4D Green’s function G 4 in eq. (20), we will now argue that 


^ ^ ((e - - 4) 


(A12) 


G2[?,e']= / d^x^G4a]J[^,e;xl], 

4e2 

where the source J obeys the equation(s) 

+ Vi + D[e]) + Vxln[fi] • Vx ((C - = 0, (A13) 

5^2 + Vi + c/[^']) J[e, C'; 4] + v± in[fi] • Vx ((e' - J[c, e'; 4 ]) = 0 , 

and the boundary condition at the spatial origin of the 21 }-x±_ plane 


= 4 ] = 1 




(A14) 


Notice, if this prescription is valid, we have attributed to the source J all the 2D Poincare symmetry 
breaking due to the presence of the potential [/; for otherwise, G 4 is itself Poincare invariant. Also, 
for some fixed elapsed time r] — rj' and fixed source spatial location x', the signal at x is going to come 
from integrating over the circle on the x_L-plane defined by \r] — r/'p = {x — x'Y' + or equivalently 
= 2a] and since it is integrated over anyway, there is no need to allow for the source J to have an 


azimuthal dependence. Note that eq. (A 12 ) provides a second means to deduce that G 2 is pure tail. 


Recalling eq. (20), 


G2[e,4 = 


dr^ 


Jo 2 


4(^-4'-4] J[^,e';4] 


(A15) 




^^We begin the derivation of equations (A13) and (A14) by highlighting that G 4 [cj] itself depends on 
spacetime solely through the Minkowski world function a. Now, apply the 2D wave operator on both 


sides of eq. (A12), interchange the order of integration and differentiation, and add and subtract V^. 

+ u) G2[^, 4 = J - Vi) G 4 


(AI 6 ) 


+ (G4 • 2d^,ad>^u + Ga ■ + Vi + [/) j) }. 


The integral in eq. (A15 I is a special case of dz5[z — a]/[ 2 :] = f[a]Q[z 2 — a]0[a — zi]; the top hat 0[z2 — a]0[a — zi] 
is the constraint that, the integral returns f[a] if a lies within the interval a G [zi, Z 2 ] and zero otherwise. 
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We have integrated-by-parts one of the onto the J, and also used the fact that G 4 depends on 
spacetime solely through a. In fact, let us observe that 


-fu • vuG4 = {-ysi • g: 


4 ) 


(A17) 


and thus the prime on the G 4 can be expressed in terms of Vj^ and then integrated-by-parts: 


• 2dfj,^ad^‘^J = 2 / d^x_LG4V_L • ( ) . 


X_L 




(A18) 


Because V_l lnx]_ = 2 x_l/x]_, 

2Vu • = V± ln[xi]a^2d<9^2 . j ^ J. 

The Green’s function equation for the 2D Laplacian 

V5_ln[x5_] = 47r(^^^^[xj_], 

together with the 4D Minkowski wave equation d^Gi = 26^^\ then allow us to gather 


(A19) 


(A20) 


+ U) G2[C, e'] = 2J[C = e-, 6±]5^^'> [? - ^'] + 47rG4[d] • C'; Ou] (A 21 ) 

+ j^ d^xxG^ { (5^20^2 + Vi + [/) J + Vu ln[xl]df,,a ■ . 

To obtain 25^‘^'> — ^'] on the right hand side, we demand that eq. (|A13|) be satisfied to have the 


second line of eq. (A21) vanish, and the simultaneous boundary conditions J[^ = C^0 _l] = 1 and — 
()_l] = 0 for all ^. But these latter two conditions amount to demanding J[^, 0_i_] = 1 


for any y, i.e., eq. (A14). 


2D Minkowski Massive Scalar We now show how this 3D picture offers an alternate perspective 


on eq. (A 8 ), the 2D massive scalar flat spacetime Green’s function G 2 , which through eq. (14), would 
yield the latter’s 4D counterpart. We have U = m?, and will further make an inspired guess, based on 


the correspondence between 2(t and ri (recall eq. (A15)), that J is actually time independent - since 


m 


2D, Poincare invariance indicates G 2 depends on spacetime solely through a. (In this specific context, 
there is a relationship between rotational O 2 invariance in the 3D picture and Poincare symmetry in 


2D.) The problem of massive waves in 2D Minkowski is now translated into, via eq. (A13), a static 
Laplace-Helmholtz type equation in 2 spatial dimensions: 


V5_ -I- ) J[r±] = 0. 


(A22) 
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(This reduces to the same equation in eq. (A9).) The solution obeying the boundary condition in eq. 
J[r_i_ = 0] = 1 is the Bessel function 




{A23) 


Inserting this into eq. (A15) yields the result in eq. (AlO). 


First derivatives on the light cone Let us employ light cone coordinates so that 

the 2D flat metric reads d.^'’'d.^“ and = Ad+d-. For some scalar field ijj obeying the homogeneous 
2D wave equation + C/)V’ = 0, suppose its value is known on the light cone based at some fixed 
location , i.e., '(/’[(^ ~ = 0] is known; this also means the first derivatives and 

along the light cone can be derived. (In fact, by this assumption, all higher derivatives 
and n > 1, can be computed.) The other first derivatives of ip 

evaluated on the light cone of is given by the following integrals: 

r^~ 

/C- 


d+i;[t = n = J f, de"“D[^+ = = ^'+,1""], 


(A24) 


9_v^[^+,r = e'-] = = e'] - 


1 


T'+ 


de"+[/[r'+,r = r]V'K"+,r = n 


(A25) 


To see this, we first observe that the relevant boundary conditions at the light cone’s apex ^ are 


satisfied. Next, we take a derivative with respect to d- on both sides of equations (A24) and 9+ on 


(A25). This recovers the wave equation d^d-tp = —(1/4)1/V’ evaluated on the light cone. 


Equations (A24) and (A25), when applied to the tail function of G 2 respecting J[a = 0] = 1 (eq. 


(A4)), hands us 


9+j[e+ = r+,n = - 


T'- 


de"-t/[^+ = r,r-i 


9_j[e+,r = n = = ri- 


(A26) 

(A27) 


The first derivatives of J at the null cone’s apex ^ is zero because J is a constant on the entire 
light cone. 

Born Series Solution for and J We now turn to an infinite Born series solution for the 
retarded Green’s function in 2D, in terms of its counterpart G^ (which obeys d^G^ = and 
the potential U. This will also yield a corresponding series for J itself. We first define the operator Q 
acting on a bi-scalar as 


(QS}l^,a^ / d^eG2[C-e]U[nsr,a (cartesian). 


(A28) 
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Then the infinite Born series is 


OO 


(A29) 


t=o 


The zeroth term is defined to be ~ ^ > 1; and written in light cone 

coordinates = rj ± x and d^^ —)• (l/2)d^'''d^“, 


c'of) i{.«'i=I (n 


'R(ia) 


d^dC 


(A30) 


X - C1-2 ]... G2+[6 - ei]C/[ei]G+[ei - 


These £-nested spacetime integrals can be further broken down into, 


Q'Gt) K,e'] = 


¥ 




(A31) 


where = 1; using the shorthand d^"’*' d^" , and for f" > 1, 

h' h' h' h' 


= [ d^CfU[^e]^e-i[^e,a 
h' 


(A32) 

(A33) 


The Born series for J itself, whenever ^ lies within the future light cone of ~ which we will denote as 

e > ?' - is 


°° (_Y 


£=0 


(A34) 


For the purposes of solving the 2D Green’s function G 2 we do not need J outside the light cone. 
Therefore, once is known for ^ lying within the future light cone of by symmetry this means 

we have obtained J[^', for inside the past light cone of by swapping the labels ^ -H- this means 
we also know for ^ lying in the past light cone of . 

Derivation Let us now justify equations (A29) through ( |A34 ). We begin by applying cl| and 
5|, to eq. (A30). For £ = £), the result is just the (5-function (5i^i[^ — For i = 1, the result is 


d‘^QGt = UGo. 


(A35) 


For £ > 1, the leftmost or rightmost integral will collapse, returning 


d'^Q^G^ = UQ^-^G^. 


(A36) 
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A direct calculation would therefore tell us 


+ 00 




(A37) 


e=o 


which is (9^ + U)G 2 = 6 , the retarded version of eq. (A1 Next, we will show that using the retarded 

flat 2D massless scalar Green’s function in this construction implies the solution for also obeys the 
retarded boundary condition. Now, for any object U, the integral 


(qg+) [e,e'] 




(A38) 


1 , 1 ) 2 2 

is zero if ^ lies outside the future light cone of , due to the presence of the 2 retarded Green’s functions 
in the integrand. We have also used the result 


G-lK - e'l = ieh - = ^e+K. f'l- 


(A39) 


When ^ the integral in eq. (A38) is really that of over the rectangular region on the 

^"-2D plane bounded by the past light cone of ^ and the future light cone of i.e., (l/4)Xi[,^,^']. It is 
zero otherwise. What we have argued, therefore, is that 




(A40) 


Let us suppose that eq. (|A31|) holds for some Then 


(q'+'G+) [^,e'] = ^ y’^^d2c±,G+[e-^,+i]t/[e£+i]X,[?£+!,-e']. (A41) 


The causality arguments employed to analyze eq. (A38) tell us, because of the two retarded Green’s 
functions, the must again be proportional to G^[^ — |']. Assuming ^ > C the integral in eq. 


(A41) becomes that of [/[^f+i]X£[^£+i, ,^']/4 over the rectangular region bounded by the past and future 


null cones of ^ and respectively; eq. (A41) is thus 


(q^+ 1G+) [^,e'] = 


(A42) 


As an aside, we remark here that the derivation up to this point is not specific to 2D. This cosmological 2D reduced 
Green’s function problem is, in fact, not the only curved spacetime example that admits a solution as an inhnite Born 
series in G, its corresponding Minkowski Green’s function, because of the “flat □’’-plus-potential form of the wave 
operator. The other instance we are aware of, is that of the minimally coupled massless scalar and photon retarded 
Green’s functions in black hole (BH) geometries put in a Kerr-Schild form [27]. There, for both Schwarzschild and Kerr 
spacetimes, the Ah term in the summation goes as (BH mass)^ times Anested spacetime integrals that are independent 
of the BH mass but involve for the Kerr case the BH’s angular momentum. 
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Recalling the integral recursion relation between Xi and in eq. (A33) then proves eq. (A31) for 
arbitrary l>\. 


Because we have shown that every term in the Born series of in eq. (A29) is proportional to 


7=y+ 


G 2 , we may in fact define J to be the latter’s coefficient, namely 

_ _ (—V 


(A43) 


i=0 


(The hrst equality is, of course, consistent with eq. (A2).) The result in eq. (A34) follows once we take 


^ to lie within the future light cone of for there G 2 [C ~ (1/2)- 


A different means of arriving at eq. (A34) is to integrate with respect to both the homogeneous 


equation = — (1/4)C/J from the light cone to some arbitrary point ^ lying within the future 

light cone of taking into account the boundary conditions J[a = 0] = 1 and d±J[^'^ = = 0. The 

result is the integral equation 

1 /■« 




(A44) 




Iterating eq. (A44) infinite times yields eq. (A34). 


2D Minkowski Massive Scalar Revisited For the massive scalar in 2D Minkowski, U = is a 


constant and every term of eq. (A34) can be evaluated. The £th term is, in light cone coordinates, 

rii 


r«3 


r?2 




(A45) 


The rightmost integral gives (^2 )• d’de one after that gives us )^/2^, 

and so on. After all i integrals are done we have ~ C~Y/■ This leads us to 

u. .n ^ rn^Y^^-e+Yir-e-Y 
to ~ ■ 

This is nothing but the infinite series representation of 


Appendix B: Why d‘^Q[a] = 4^(2) 


The primary goal in this section is to justify eq. (A6). Via a direct calculation, with = (^ — Y)r- 
and — ^'Y = 2d', 


cl^0[(j] = 2da = 0. 


(Bl) 
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Therefore 0 [(t] is the homogeneous solution to the 2D massless scalar wave equation almost everywhere. 
Let us employ light cone coordinates, 


>~l~ _ I f I f 

4 — r] ±x, 4 — r] ±x , 


(B2) 


so that a = (l/2)(^+ — — i'~) and 9^ = Ad+d-. Note that ©[ct] is a constant inside the light 

cone, so d'^Q[a] can only be non-zero on the light cone, when either or . To argue 

that it is in fact only non-trivial at the apex of the light cone, i.e., ^ we shall first let be 

arbitrary, but integrate the left-hand-side of eq. (A6) about 
/■?'"+o+ 


dr45_a+0[a] = (5[e+ - C'+]4 • sgn[r - 


/g'--0+ 


= 85[e+ - e 


'+1 


(B3) 


We have used the distributional identity (5[a6] = 5[6]/|a|. Similarly, 

/■5'++0+ 


/^'+-o+ 


d^+4a+5_0[n] = 85[r - ? ]• 


(B4) 


Therefore, (9^0[a] is non-zero only when both —?■ and —?■ simultaneously. This justifies 

the (written in Cartesian coordinates) on the right hand side of eq. (A6). The factor of 4 

can now be checked by re-expressing the (5-functions of eq. (A6) in light cone coordinates, (5^^^[^ —= 
2(5[^'*“ — C^'’~](5[C~ ~ then integrating it with respect to d^'’'d^“ over an infinitesimal region 


around ^ to obtain 8. On the other hand, eq. (B3) is already the 5^0 integrated with respect 


to d^ , so to recover this 8 one simply has to integrate over d^"*". The same check can be made by 


integrating both sides of eq. (B4) with respect to d^ 
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